CLNS 00-1695 
CLEO 00-19 



Observation of the Charmed Baryon at CLEO 



CLEO Collaboration 
(February 7, 2008) 

Abstract 



The CLEO experiment at the CESR collider has used 13.7 fb _1 of data 
to search for the production of the (ess-ground state) in e + e~ colli- 
sions at y/s ~ 10.6 GeV. The modes used to study the are Q~ir + , 
il~7T + 7r°, "E~K~it + it + , H°A'~7r + , and f2~7r + 7r + 7r~. We observe a signal of 
40.4±9.0(stat) events at a mass of 2694.6±2.6(stat)±1.9(syst) MeV/c 2 , for 
all modes combined. 
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Many experimental groups have searched for the in numerous decay modes: how- 
ever, their reported masses are only marginally consistent with each other. The WA62 
experiment JI[ claimed the first evidence of in the H~ K~k + -k + decay mode with a mass 
of 2740.0 ± 20.0 MeV/c 2 . The ARGUS Collaboration @ published an signal in the 
Z-K-7T+7V+ mode with a mass of 2719.0 ± 7.0 ± 2.5 MeV/c 2 , based on 0.380 fb" 1 of data. 
This result was contradicted by CLEO, in an unpublished conference paper |J using 1.8 
fb- 1 of data. Later, E687 |J published an tt° c mass of 2705.9 ± 3.3 ± 2.3 MeV/c 2 using the 
f2~7r + mode and a mass of 2699. 9± 1.5±2.5 MeV/c 2 using the higher-statistics Y, + K~K~7T + 
mode. In 1995, the WA89 Collaboration |5j reported 200 events in seven decay modes 
with an average mass of 2707. Oil. 0(stat) MeV/c 2 ; this result remains unpublished. 

The (c{ss}) is a J p = | + ground state baryon, where {ss} denotes the symmetric 
nature of its wave function with respect to the interchange of light-quark spins. Various 
theoretical models [§-|9| predict an mass in the range 2664 - 2786 MeV/c 2 . 

The data used in this analysis were collected with CLEO II |I(J and the upgraded 
CLEO II. V |y] detector operating at the Cornell Electron Storage Ring (CESR). The data 
corresponds to an integrated luminosity of 13.7 fb _1 from the T(4S) resonance and at 
energies in the continuum region just below. We searched for the in the five decay modes 
Q~7i + , f2~7T + 7r°, Q~7T + 7T + 7T~ , H~fT~7r + 7r + , and E°K~"tt + . The choice of these five modes is 
based mainly on the pattern of other charmed baryon decays. Reconstruction efficiencies 
and the size of the combinatorial background were other considerations. A sixth channel, 
T l + K~K~n + , was also investigated because E687 [Q showed a significant signal in this decay 
mode, although CLEO has rather low efficiency for this mode. 

Charmed baryons at CESR are either produced from the secondary decays of B mesons 
or directly from e + e~ annihilations to cc jets. We introduce x p as the scaled momentum 



of a candidate, where x p = p/p ma x, and p max = y E^ — M 2 with E^ equal to the beam 
energy and M the mass of the candidate. Our search is limited to x p > 0.5 or x p > 0.6, 
depending on the decay mode, so as to avoid the combinatorial background that dominates 
at low x p . Charmed baryons from B meson decays are kinematically limited to x p < 0.4, so 
our search is limited to the Vl° c baryons produced in the e + e~ continuum. We implemented 
p/K/tt identification by defining a joint probability for each hypothesis, using both the 
specific ionization (dE/dx) in the drift chamber and the time-of-flight to the scintillation 
counters. A charged track is defined to be consistent with a particular particle hypothesis 
if the corresponding probability is greater than 0.1%. We required all the charge tracks in 
all the decay modes to be consistent with their respective particle hypotheses. To further 
reduce the combinatorial background we also required the momentum of daughter pions and 
kaons from f2° to be greater than 0.2 GeV/c to 0.5 GeV/c depending on the decay mode. 

We begin the analysis by reconstructing A— >-p7r _ , H°^A7r°, H~^A7r~, Q~^AK~ , and 
E + ^p7r°. Charge conjugation is implied throughout the analysis. The analysis procedure for 
reconstructing these particles closely follows that presented elsewhere [12-14]. The hyperons 
are required to have vertices well separated from the beamspot, with the flight distance of 
the secondary A greater than that of the H°, or Q~ . We then combine these hyperons 
with tracks from the primary event vertex to reconstruct candidates. 

In each mode the signal area above the background is obtained by fitting with the sum 
of a Gaussian signal function (with width fixed at the signal Monte Carlo predicted value 
for that mode) and a second order polynomial background. The Monte Carlo sample used 



4 



in this analysis was generated for the two CLEO detector configurations using a GEANT- 
based jn| simulation and was processed similarly to the data. We simultaneously fit the 
five modes to a single mean value for the mass. In the Q~tt + mode, we required x p to be 
greater than 0.5 and the 7r + momentum to be greater than 0.5 GeV/c. Figure |l|(a) shows 
the invariant mass distribution; a fit to this distribution yields a signal of 13.3 ±4.1 events. 
In the f2~7r + 7r° mode, only 77 combinations having invariant mass within 12.5 MeV/c 2 (2.5 
standard deviation) of the nominal 7i~° mass are used as tt° candidates; we assume the 
photons used for reconstructing 7r°^77 come from the event vertex. We also required x p 
to be greater than 0.5 and the tt + and tt° momenta to be greater than 0.3 and 0.5 GeV/c, 
respectively. Figure |I](b) shows the invariant mass distribution. A fit to the distribution 
gives a signal yield of 11.8 ± 4.9 events. Figure |l|(c) shows the Q~7t + 7t~tt + invariant mass 
distribution for x p greater than 0.5. All the charged pions are required to have momenta 
greater than 0.2 GeV/c. The corresponding fit yields a signal of —0.9 ± 1.4 events. In the 
E°K^tt + mode, we considered combinations with x p greater than 0.6, since the combinatorial 
background is higher in this mode. Figure |I|(d) shows the invariant mass distribution, with a 
fit yielding a signal of 9.2±4.9 events. In the H~ K~7t + 7t + mode, we required x p to be greater 
than 0.6 and pion and kaon momenta to be greater than 0.2 and 0.3 GeV/c, respectively. 
A fit to this distribution yields a signal of 7.0 ± 3.7 events. Finally, in the Yi + K~ K~tt + 
mode, we required x p to be greater than 0.5 and required charged track momenta to be 
greater than 0.3 GeV/c. The upper limit on the yield is 9.5 events (90 % C.L.). Figure |l|(f) 
shows the invariant mass distribution for Y< + K~K~tt + mode. The efficiency for Yi + K~ K~tt + 
reconstruction is ~ 15% of that for the Q~tt + mode, which has the highest signal yield. We 
have not included the E + K~ K~7t + mode in the mass measurement. The total yield in the 
five decays modes, excluding I] + K~K~7t + , is 40.4 ± 9.0 events as shown in Table |. The 
corresponding combined mass distribution is shown in Figure ^|. 

To better determine the mass, we have performed an unbinned maximum-likelihood 
fit using the sum of a single Gaussian and a second order polynomial background. There 
are two inputs to the fit, the invariant mass Mi and the corresponding mass resolution Oi 
of each of the 458 mass candidates from 2.55 to 2.85 GeV/c 2 . The invariant mass Mj of 
each candidate is calculated using a vertex constrained fit; the uncertainty <7j is obtained 
from the covariance matrix of the fit. The likelihood function to maximize is the product 
of probability density functions (PDFs) for all the candidate events, and has the following 
form: 

£(M(fi°),/ 8J ai,a 2 ) = n 

i 

where G(y\a) = (l/V2^)exp(-y 2 /2a 2 ) and P{y) = 1.0 + a x {y - 2.7) + a 2 {y - 2.7) 2 . 
M(f2°) is the fitted mass, S is a global scale factor multiplying cr^, and f s is the fraction of 
signal events. We tested the fitting procedure by first applying it to Monte Carlo generated 
events with the mass set to 2695 MeV/c 2 . The fitted mass using the above PDF is 
2694.9±0.1(stat) MeV/c 2 for the Monte Carlo and 2694.6±2.6(stat) MeV/c 2 for the data. 
The fitted scale factor 5" is 1.72±0.42 for the data and 1.16±0.02 for the simulated events. 
The fitted value of f s = 0.099±0.027(stat) implies a signal size of 45.2±13.7(stat) in good 
agreement with the result from Table |. We tested our fitter on the S° charmed baryon 
in the S _ 7r + final state. Our fit gives a 5° mass that is consistent with the world average 



f s G(Mi-M(n° c )\Sai) + (l-f 8 
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value, and the fitted global scale factor S for the S° Monte Carlo and the data are in good 
agreement with each other. 

We have also checked for goodness-of-fit by performing a binned-likelihood version of the 
above fit (but without event-by-event mass uncertainties). This gives an almost identical 
mass value and similar yields, and forms the basis for the curve shown in Figure 0. The \ 2 
for the fit to the combined data is 46.2 for 46 degrees of freedom. 

The dominant systematic uncertainty in the mass measurement comes from its sensitivity 
to the fitting method employed - the difference in the weighted average of the individually 
fitted modes and the unbinned maximum likelihood method. The various fitting methods 
contribute 1.5 MeV/c 2 to the systematic uncertainty. A systematic uncertainty (0.5 MeV/c 2 ) 
comes from imperfect treatment of modes with tt° mesons due to mismeasured photons at low 
energies that give rise to an asymmetric ir° peak. Additional contributions to the systematic 
uncertainty come from uncertainties in the energy loss correction for charged tracks (0.25 
MeV/c 2 ) and the overall momentum scale (1.0 MeV/c 2 ). Taking these errors in quadrature, 
we estimate a total systematic uncertainty of 1.9 MeV/c 2 . 

In Table | we also give the measured inclusive cross-section times the branching faction, 
a-B, for x p > 0.5 into ir7r + , fi-vr+Tr , E-K-n + 7r + , E°K'tt+, fi-7r+7r+7r~and T I + K~K-'K + . 
We estimated the systematic errors for the branching fractions by changing the mass by 
±1.0cr (combined error) from its best fit value. In Table | the first error is due to statis- 
tics and the second, when given, to systematics. In charm decays, c— >Ws, the W tends 
to couple more strongly to two pions (via the p-meson) than to a single pion. The rela- 
tive branching fraction B(c— >s7t + tt°)/B(c— >stt + ) is greater than unity. The relative branch- 
ing fractions for B(D + ^Wn + n°) /B(D + ^Wn + ), B(D ^K-tt + tt°)/B(D°^K-tt + ), and 
£(A+^A7r+7r°)/i3(A+^A7r + ) are 3.4±l.l(stat), 3.6±0.2(stat), and 4.0±1.9(stat) |]|, re- 
spectively. We also observe a similar trend in the branching fraction for fl®— >Q~7t + tt° rel- 
ative to f2°-^f2~7r + , as given in Table |I[ We have also studied the momentum spectrum of 

finding it consistent with other charmed baryons WJft. 



TABLE I. fig results in different decay modes. The fourth column shows the branching 
fractions relative to f2~7r + ; the fifth column shows the cross section times branching fraction. The 
fourth column has only statistical, while fourth and fifth columns have statistical and systematic 
uncertainties. 





a MC (MeV/<?) 


Fitted Yield 


Relative B 


a ■ B (fb) 






mode dependent x p 


all x p > 0.5 


all x p > 0.5 




5.87 


13.3±4.1 


1.0 


11.3 ±3.9 ±2.0 


tt-TT+TT 


9.71 


11.8±4.9 


4.2 ±2.2 ±0.9 


47.6 ± 18.0 ±3.1 


E°K-tt+ 


6.72 


9.2±4.9 


4.0 ±2.5 ±0.4 


45.1 ±23.2 ±3.7 


H~ K~ 'ir^ ~ir + 


5.46 


7.0±3.7 


1.6 ± 1.1 ±0.4 


18.2 ±10.6 ±3.3 




4.89 


-0.9±1.4 


< 0.56 


< 5.1 @ 90 % CL 


Combined 5 modes 




40.4±9.0 






Y,+K-K-it + 


6.18 


2.8±4.1 


< 4.8 


< 53.8 @ 90 % CL 



In conclusion, we observe the with a mass of 2694.6±2.6(stat)±1.9(syst) MeV/c 2 in 



6 



the five decay modes f2~7r + , f2~7r + 7r°, f2~7r + 7r + 7r~, S~ K~tt + tt + , and E K~tt + . Although 
the signal is not statistically significant in any individual mode, the combined signal stands 
out over the background with a yield of 40.4±9.0(stat) events. Our measured a ■ B value for 
the E~K~n + ir + mode (18.2 ± 10.6 ± 3.3 fb) is in clear disagreement with the ARGUS value 
(2410 ± 900 ± 300 fb) for the same x p range. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, the Natural Sciences 
and Engineering Research Council of Canada, the A. P. Sloan Foundation, the Swiss Na- 
tional Science Foundation, the Texas Advanced Research Program, and the Alexander von 
Humboldt Stiftung. 
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FIG. 1. The mass distribution and simultaneous fit to the five modes: (a) Q 7r + , (b) 
irvr+Tr , (c) n-ir+ir+ir-, (d) Z°K-ir + , (e) ~.-K-it + it + . The mode (f) T i + K-R-tt + has not 
been included in the fit. The signal is fit with a gaussian of fixed width while the background is 
fit to a second order polynomial. 
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FIG. 2. The invariant mass distribution for the sum of f2 ir + , f2 ir + ir , f2 ir + ir + ir ,S°K ir 
and H _ iT _ 7r + 7r + conibinations. The fit function is a sum of the fit functions from Figure |l|. 
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